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Abstract 

The aim of this paper is to classify the bispectral operators of any rank 
with regular singular points (the infinite point is the most important one). 
We characterise them in several ways. Probably the most important result 
is that they are all Darboux transformations of powers of generalised Bessel 
operators (in the terminology of |BHY1]). For this reason they can be effectively 
parametrised by the points of a certain (infinite) family of algebraic manifolds 
as pointed out in [BHY1]. 
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Introduction 

The present paper is devoted to the characterisation and the classification of bispec- 
tral operators of any rank and order with only regular singularities. Before stating 
our results and placing them properly amongst the other research we would like to 
give few definitions and to recall some of the fundamental results in the area. 

An ordinary differential operator L(x, d x ) is called bispectral if it has an eigen- 
function i/j(x,z), depending also on the spectral parameter z, which is at the same 
time an eigenfunction of another differential operator A(z,d z ) now in the spectral 
parameter z. In other words we look for operators L , A and a function ifi(x,z) 
satisfying equations of the form: 



14 = 

Aift = 9(x)tp. 



(0.1) 
(0.2) 



Initially the study of bispectral operators has been stimulated by certain problems 
of computer tomography (cf. [Gl, G2]). Later it turned out that the bispectral 
operators are connected to several actively developing areas of mathematics and 
physics - the KP-hierarchy, infinite-dimensional Lie algebras and their representa- 
tions, particle systems, automorphisms of algebras of differential operators, etc. (see 
e.g. |PHY1|, PHY4 |BHY5|, |BW|, pGl, W% |W|, as well as the papers in the 
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proceedings volume of the conference in Montreal [|BP|| ). There are also indications 



for eventual connections with non-commutative algebraic geometry [W3 
In the fundamental paper 



Duistermaat and Griinbaum raised the problem 
to find all bispectral operators and completely solved it for operators L of order two. 
The complete list is as follows. If we present L as a Schrodinger operator 



L 



the bispectral operators, apart from the obvious Airy (u(x) = ax) and Bessel (u(x) = 
cx~ 2 ) ones, are organised into two families of potentials u(x), which can be obtained 
by finitely many "rational Darboux transformations" 

(1) from u{x) = 0, 

(2) from u{x) = — {\)x~ 2 . 



Thus the classification scheme prompted by the paper [DG] is by the order of 



the operators. G. Wilson [Wl] introduced another classification scheme - by the 
rank of the bispectral operator L. We recall that the rank of the operator L is 
the dimension of the space of the joint eigenfunctions of all operators commuting 
with L. For example all the operators of the family (1) have rank 1, while those of 
the family (2) have rank 2. 



In the above cited paper |Wlfl (see also |W2| ) Wilson 
gave a complete description of all bispectral operators of rank 1 (and any order). In 
the terminology of Darboux transformations (see | BHY1| 1) all bispectral operators of 
rank 1 are those obtained by rational Darboux transformations on the operators with 
constant coefficients, i. e. polynomials p(d x ). Several beautiful connections of the 
bispectral operators with KdV and KP-hier archies, algebraic curves and Calogero- 



Moser particle systems have also been revealed in [ DG , Wl 



We will not touch upon all results in the papers [DG, WL] but we would like to 
point that in both of them the classification is split into two, more or less independent 
parts. First, there is an explicit construction of families of bispectral operators of a 
given class (order 2 in [ DG| ; rank 1 in | Wlj ) The construction can be given in terms 
of Darboux transformations of "canonical" operators. The second part is to give a 
proof that, if an operator (in the corresponding class) is a bispectral one, then it 
belongs to the constructed families. 

In several other papers devoted to the bispectral problem (see [0, [G^, KRo| ]) 
the authors deal with an analog of the first part of the problem, i.e. they construct 
new families of bispectral operators. The most complete results in that direction 



have been obtained in [BHY1, BHY4 |. To the best of our knowledge, all known up 
to now families of bispectral operators can be constructed by the methods of the 
latter papers. A challenging problem is to prove that all the bispectral operators 
have already been found. A natural approach would be to divide the differential 
operators into suitable classes , e. g. - by order as in DG ] or by rank and to try to 
isolate the bispectral ones amongst them. But having in mind the constructions of 



the fundamental papers [DG, Wl], with their different and quite involved methods, 
the complete classification seems to be a difficult and lengthy project. One may try 
to consider the operators with a fixed type of singularity at infinity. Obviously, then 
there arises another difficult problem - to determine what restrictions on the kinds 
of singularities are imposed by the condition of bispectrality. 



2 



In the present paper we consider the class of operators with regular singularities 
at infinity. In fact the main results sound much stronger. To explain them we 
introduce some definitions and notations which will be used also throughout the 
paper. We are going to consider operators, normalized as follows: 



N 



L = Y J Vk(x)d k x , (0.3) 

where the coefficient at the highest derivative Vjv = 1 an d the next coefficient 
V/v-i = 0. Now our assumption is that 

limV 7 -(x)=0, j = 0, ...,N — 1 when x — » oo. (0-4) 

(It is well known that with the above normalization all coefficients of L are rational 
functions (see DG, Wl|) and hence (|0.4|) makes sense). 



Important examples of such operators are the generalized Bessel operators. As 
we are going to use them throughout the paper we recall the definition. Introduce 
the notation D = xd x . 

Definition 0.1 Generalized Bessel operators Lg are the operators 

L f3 = x- N (D-(3 1 )...(D-(3 N ), (A,..,fe)6C w , (0.5) 

In what follows we will call the above operators by abuse of terminology (but for 
simplicity) Bessel operators. 

After this preparation we can formulate the result which is the core of the present 
paper. 



Theorem 0.2 Let L be a bispectral operator ( p.3[ ) with coefficients satisfying ( p.4[ ). 
Then L is a monomial Darboux transformation of a Bessel operator. 

The class defined by ( |0.3j ) and (p.4| ) includes essentially all the bispectral oper- 
ators found in [DG|: the Bessel operators and both of the families (1) and (2), the 



only exception being the Airy operator. On the other hand it includes one of the 



most interesting classes , found in [BHY1|. These are the operators obtained by 



Darboux transformations on powers of the Bessel operators. This class was later 
characterized as follows. In [ BHY2| ] there have been constructed highest weight 



modules Aig with highest weight vectors - the corresponding to ( p.5| ) r-functions 



Tg. Then in [BHY5| it is shown that the r-functions in the modules A4g are exactly 



the T-functions of the operators which are monomial Darboux transformations. 



In the course of performing the proof of Theorem 0.2 we show that the assump- 
tions Q0.3j ) and ( p.4[ ) for the bispectral operator L impose further restrictions on it, 
which justify partially the title. 

Theorem 0.3 If the bispectral operator ( |0,3j ) satisfies ( p.4| ), then the point x = oo 
is a regular singular point. 
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The proof of this theorem is probably the most involved part of our constructions 
(see Sect. 3). The regularity of the finite points follows indirectly from Theorem 



In Sect. 4 we give another characterization of the bispectral operators ( |0.3[ ) with 
the restriction ( p.4[ ). 



Theorem 0.4 Any rank r bispectral operator L is Z r -invariant. 



The result is interesting and natural by itself (cf. [ pG| , BHY1| ]) but in the present 
paper it is also the next step in our final goal. 

Finally in Sect. 6, putting together the different pieces of our construction in 



the preceding sections and using the main results of [BHY1, BHY5] we obtain the 



following complete characterization of the Fuchsian bispectral operators. 

Theorem 0.5 The following conditions on the operator L in the form ( |0.3| ) are 
equivalent: 

1) L is bispectral and satisfies ( |0.^j ); 

2) L is bispectral and has only regular singular points (i.e., L is Fuchsian); 

3) L is a monomial Darboux transformation of a Bessel operator ( |0.-5| ); 

4) The corresponding to L t- function belongs to one of the modules Mp. 

In the case when the order of L is two the equivalence between 1) and 3) contains two 



of the most important (and difficult) theorems of [DG], concerning the families (1 



and (2) above. In that sense the present paper represents their direct generalization. 



The methods which we utilize have some resemblance to the ones used in [DG]. 
In particular the Darboux transformations constitute one of the main steps of our 
proof. But as a whole we use different ideas. First, we work essentially with the 
algebraic structure of different rings of differential or pseudo-differential operators. 
Essentially we do not use the wave function as in |DG| ]. This we achieve by using 
the bispectral involutions on pseudo differential operators in Sect. 2. In the same 
section we observe that a bispectral operator L (with the restrictions fl0.3| ) and ( p.4| )) 
satisfies a variant of the so called "string equation" : 

[L,Q}=NL n+ \ (0.6) 



where Q is an operator built out of L. The equation fl0.6|) prepares us to use 
certain techniques from differential algebra in order to study the singular point 



of L at infinity. In particular we use the methods invented by J. Dixmier [Dx] 
in his studies on the Weyl algebra. Roughly speaking one associates with each 
differential operator L a quasi-homogeneous polynomial pi(X,Y) in such a way 
that it contains the information about the "worst" terms of L (in our case these 
are the most irregular ones). See [|Dx| and Sect. 3 for more details. Then in the 
same section the analysis of pi(X, Y) shows that the assumption of irregularity of 
the point x = oo is incompatible with the string equation ( |0.6| ). 

The techniques from Sect. 2 is used also in Sect. 4 to prove that the rank r of 
the operator L imposes its Z r -invariantness. Using it and the fact that the infinite 
point is regular it is easy to perform Z r - invariant Darboux transformations on L in 
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order to reduce the number n in the string equation (0.6) to 0. This automatically 



gives that the operator obtained in this way is a Bessel operator. 

At the end of the introduction we point out that our method treats all ranks 
and orders in one scheme. We expect that some of its components can be useful in 
other classification problems. 

Aknowledgements. We are grateful to L. Haine and G. Wilson for a number 
of discussions and suggestions that lead to considerable improvements of the paper. 
The first mentioned author acknowledges the kind hospitality of Universite de Lou- 
vain, Universite de Nantes and Universite Toulouse III, where the last part of the 
paper was done. 



1 Preliminaries 

In this section we have collected some terminology, notations and results relevant for 
the study of bispectral operators. Our main concern is to introduce unique notation 
which will be used throughout the paper and to make the paper self contained. 
There are also few results which cannot be found formally elsewhere, but in fact are 
reformulations (in a suitable for the present paper form) of statements from other 
sources. 



1.1 

In this subsection we recall some definitions, facts and notation from Sato's theory 
of KP-hierarchy H DJKM, 3W| needed in the paper. For a complete presentation 
of the theory we recommend also fDi| vM|. We start with the notion of the wave 
operator K(x,d x ). This is a pseudo-differential operator 

oo 

K(x,d x ) = l + jraj{x)&, (1.1) 
i=i 

with coefficients dj(x) which could be convergent or formal power (Laurent) series. 
In the present paper we will consider a,j most often as formal Laurent series in 
The wave operator defines the (stationary) Baker- Akhiezer function i{j(x,z): 

*P(x,z)=K(x,d x )e xz . (1.2) 

^From ([O]) and ( |l.2[) it follows that tp has the following asymptotic expansion: 

oo 

ip(x,z) = e xz (l + Y, a i(x) z ~ j ), z^oo. (1.3) 
l 

Introduce also the pseudo-differential operator P: 

P(x,d x ) = Kd x K- 1 . (1.4) 

The following spectral property of P, crucial in the theory of KP-hierarchy, is also 
very important for the bispectral problem: 

Ptp(x, z) = zip(x, z). (1-5) 
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When it happens that some power of P, say P , is a differential operator, we get 
that tp(x, z) is an eigenfunction of an ordinary differential operator L = P N : 



Lip = z N ip. 



(1-6) 



It is possible to introduce the above objects in many different ways, starting with 
any of them (and with other, not introduced above). For us it would be important 
also to start with given differential operator L: 

L(x, d x ) = + V N . 2 (x)d N - 2 + . . . + V (x). (1.7) 

One can define the pseudofferental operator P as an A^-th root of the operator L: 

P = L7T=d+..., (1.8) 

and the wave operator K as: 



LK = Kd 



N 



[1.9) 



An important notion, connected to an operator L is the algebra Al of operators 
commuting with L (sec | |Kr| , PQ |). This algebra is commutative one. The wave 
function tp(x,z) (defined in (|1.2j)) is a common wave function for all operators M 
from Al' 

MiP(x,z)=g M (z)il>(x,z). (1.10) 

We define also the algebra Al of all functions gu{z) for which ( l.lOj ) holds for some 
M £ Al- Obviuosly the algebras Al and Al are isomorphic. 



Following [Ki] we introduce the rank of the algebra Al as the greatest common 
divisor of the orders of the operators in Al- 

1.2 

Here we shall briefly recall the definition of Bessel wave function and of monomial 
Darboux transformations from it. For more details see [BHY1]. Let f3 6 be such 
that 

X N(N — 1) 



i=i 



(i.ii) 



Definition 1.1 [pj |2| |BHYl|1 Bessel wave function is called the unique wave func- 
tion ^p(x, z) depending only on xz and satisfying 



L p (x, d x )^ p {x, z) = z N ^ p {x, z), 

where the Bessel operator Lp(x,d x ) is given by (p^q). 

Because the Bessel wave function depends only on xz, ( |1,12| ) implies 

D x V[)(x, z) = D z V a (x,z), 
L p {z, d z )Vp(x, z) = x N Vp(x, z). 



(1.12) 



(1.13) 
(1.14) 



To introduce the monomial Darboux transformations of Bessel wave functions we 



first recall the definition of polynomial Darboux transformations given in | BHY1 ]. 
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Definition 1.2 We say that the wave function ^ is a Darboux transformation of the 
Bessel wave function ^p{x,z) iff there exist polynomials f(z), g{z) and differential 
operators P(x,d x ), Q(x,d x ) such that 

® = -^-P(x,d x )Vp{x,z), (1.15) 
*p(x,z) = j^rQ(x,d t )*. (1.16) 

The Darboux transformation is called polynomial iff the operator P(x,d x ) from 
( |1.15| ) has the form 

P(x,^) = x- n ^p fc (x 7V ) J D x fc , (1.17) 

A:=0 

where are rational functions, p n = 1. 

We will need the following two definitions of monomial Darboux transformations. 
Their equivalence is proved in [BHY1]. 

Definition 1.3 We say that the wave function ^(x, z) is a monomial Darboux trans- 
formation of the Bessel wave function ^p(x, z) iff it is a polynomial Darboux trans- 
formation of *£p(x, z) with g(z)f(z) = z dN , d E N. Further the differential operator 

L = d M + V M -2d M - 2 + ... + V 

is a monomial Darboux transformation of Lp if the wave function corresponding to 
L is a monomial Darboux transformation of the wave function corresponding to Lp. 



Definition 1.4 The wave function ^f(x, z) is a monomial Darboux transformation 
of the Bessel wave function ^p(x,z) iff fll.15 ) holds with g(z) = z n , n = ordP and 
the kernel of the operator P(x, d x ) has a basis consisting of several groups of the 
form 

k mult(Pi+kN)-l 



it 



fc=0 



i=o 



y=ln x 



< / < jo, 



(1.18) 



where mult(/?j + kN) := multiplicity of /3$ + kN in Uj=i{/?7 + NZ>q} and jo 
max{j|6fcj 7^ for some k}. 



^From Definition L2 and Definition L3 one immediately obtains the following de- 
scription of monomial Darboux transformations: 

Lemma 1.5 The differential operator L is a monomial Darboux transformation 
of the Bessel operator Lp iff there are differential operators P = P(x,d x ) , Q = 
Q(x,d x ) and numbers d , d! such that 



Q{x, d x )P{x, d x ) = Lp{x, d x ) d , 

P(x,d x )Q(x,d x )=L(x,d x f, 
where the operator P satisfies (|1.17|) . 



(1.19) 
(1.20) 



7 



We will also reformulate some results from [BHY1]. In [BHY1| one can find a proof 
of the following statement. 



Lemma 1.6 If Lg is a Bessel operator of order N and rank r, there exists a Bessel 
operator Lg< of order r such that Lg is a monomial Darboux transformation of Lg< . 



For the proof of this lemma see the proof of Proposition 2.4 from [ BHY1| (although 



the statement there is formulated in a different way). We end this subsection by 
reformulating (in a weaker form) the main result, which we need from [BHY1]. 



Theorem 1.7 The monomial Darboux transformations of the Bessel operators are 
bispectral operators. 

1.3 



Here we recall several simple properties of bispectral operators following [DG, Wl]. 



As we have already mentioned in the introduction we are going to study ordinary 
differential operators L of arbitrary order N which are normalised as in ( p. 3D , i.e. 
with V/v = 1 and V/v-i = 0. Assuming that L is bispectral means that we have also 
another operator A, a wave function iff(x, z) and two other functions f(z) and 9(x), 



such that the equations (|0.l| ) and (P-2j ) hold. The following lemma, due to |DG[| , has 
been fundamental for all studies of bispectral operators. 

Lemma 1.8 There exists a number m, such that 

(ndL) m+1 9 = 0. (1.21) 



For its simple proof, see [DG, Wl|. We will consider that m is the minimal number 
with this property. An important corollary of the above lemma is the following 
result. 

Lemma 1.9 The functions f(z) and 9(x) are polynomials. 



The next result is also contained in jDGL Wl | , but it is not formulated as a separate 



statement. We give its short proof following Wl] 



Lemma 1.10 The coefficients ctj in the expansion ( |1.1[ ) of the wave operator K are 
rational functions. 

Proof. ^From the equation ( |1.21| ) it follows that 

(addfr+^irtK) = 0. 

On the other hand the kernel of the operator (add^) m+1 consists of all pseudo- 
differential operators whose coefficients are polynomials in x of degree at most m. 
This gives that 

6K = KQ, (1.22) 
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with a pseudo-differential operator G : 

oo 

e = e + X>;V (i-23) 
i 

whose coefficients @j are polynomials of degree at most m. We have 6 = Qq. 
Comparing the coefficients at dx we find that all the coefficients aj(x) of K are 
rational functions. □ 

Remark 1.11 We notice that at least one of the coefficients of Qj has degree exactly 
m, where m from Lemma is minimal. This fact will be used later. 

The last lemma has as an obvious consequence one of the few general results, im- 
portant in all studies of bispectral operators. Noticing that the coefficients of L are 
polynomials in the derivatives of aj(x) we get 

Lemma 1.12 The coefficients of L are rational functions. 



2 Bispectral involutions and the string equation 

The condition ( p.4| ) for vanishing of the coefficients Vj(x) of a bispectral operator L 
implies further restrictions on all objects connected to L - the wave function ip(x, z), 
the wave operator K and the coefficients of L itself. This gives us the opportunity 
to define two anti-isomorphisms b and b\ ("bispectral involutions") between the al- 
gebras of of pseudo-differential operators with coefficients - formal Laurent series 
in the variables x~ l and z . In its turn using these anti-isomorphisms will allow 
us to continue our further constructions in the rest of the paper by purely alge- 
braic analysis on the differential or pseudo-differential operators, avoiding the wave 
function. 

2.1 Bispectral involutions 

In the next lemma, following [DG| we find the simplest restrictions on the coefficients 



of the wave operator K and on L. 

Lemma 2.1 (i) The coefficients Vj(x),j = N — 2, ...0 of L vanish at oo at least as 
(ii) The coefficients ctj, j = 1, . . . of the wave operator K vanish at least as 



x 2 



x 1 . 



Proof. We are going to prove both statements simultaneously. We use the formula 

LK = Kd N , 



Lemma l.ldj and Lemma 1.12. Comparing the coefficients at d N 2 at the both sides 



of the above identity we get: 

V N -2 + Na[ = 0. 
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Having in mind that Vjv-2 is equal to the derivative of the rational function ot\ and 
that it vanishes at oo we see that it vanishes at least as x~ 2 . Continuing in the same 
manner we find 

N(N-l) n , 
V/v-3 + V N - 2 cti + — ^ -«i + Na 2 = 0. 

We see that a' 2 is vanishing (at least as x~ 2 ) and that Vn-3 vanishes again at least 
as x~ 2 , being a sum of such terms. By induction we get that o; s _ 1 , s = 1, ...N — 1 
vanishes at least as x~ 2 and the same holds for Vat-s, s = 2, .., N as it is a sum of 
products VjCtm\ where N — 1 > j > s, m = 1,...,N — 1 (here denotes k-th 
derivative) and also pure derivatives of a m . Arguing as above we get the statement 
of the lemma. □ 

Following [BHY4] we will introduce an anti-isomorphism b between the algebra 



B of pseudo-differential operators P{x,d x ) in the variable x and the algebra £>' of 
pseudo-differential operators R(z, d z ) in the variable z. More precisely B consists of 
those pseudo-differential operators 



CO 



k 

for which there is a number n G Z (depending on P) such that x n pj(x~ 1 ) 
k, k + 1, . . . are formal power series in x~ l . The involution 

b : B — >B' 



is defined by 

oo 

b(P)e xz = Pe xz = Y,z~ j P j (d; 1 )e xz , for PeB (2.1) 

k 

i. e. b is just a continuation of the standard anti-isomorphism between two copies 
of the Weyl algebra. In what follows we will use also the anti-isomorphism 

bi:B — ► B', b 1 {P) = b(Ad K P). (2.2) 

Obviously b and b\ can be considered as involutions of B and without any ambiguity 
we can denote the inverse isomorphisms b^ 1 , b^ 1 : B' — ► B by the same letters. 

Remark 2.2 If we use relations ( p.l| ) and ( p.2| ) to define an involution b\ on the 
subalgebra of B generated by L and 9, then we have 

h(L) = b(K~ l LK) = b(Ad K L), 

bx(0) = biK^QK) = b{Ad K 0). 
This prompts the definition ( |2.2[ ). 
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Since the operators K and 6 are from B we can define two operators S and A as 
follows: 



S(z,d z ) = b(K(x,d x )), (2.3) 
A(z,d z ) = 6(6). (2.4) 



Explicitely one has 



s = J>- J «;(&) = $>;(*" W, a = 1 (2.5) 

i=o i=o 

and also 

oo m 

A(z,d z ) =J2z~ j ®j(dz) = £A < (*- 1 )^, (2.6) 

j=0 i=0 

where A m 7^ (see Remark |1.11| ) and the coefficients Aj and a,j should be viewed as 
formal power series. We are going to see that they are polynomials in z . 

Lemma 2.3 The coefficients aj of the operator S are polynomials in z" 1 . 

Proof. Using that 

LK = Kd N , 
we can apply the involution b and to derive: 

Sb(L) = z N S. 
Rewrite in details the last formula: 



oc 



(Eajiz- 1 )^)^ + z N - 2 V N . 2 (d z ) + ...) = z N (Y j a J (z~ 1 )d^). 


Comparing the coefficients at dg 3 for j = 2, 3... and having in mind that according 
to Lemma 2.1: 



V k (d z ) = J2V k;S d- s , k = 0,...,N-2, 



2 

we obtain relations for a\ and a 2 in the form: 

N-2 



-Nz N ~ l ai + J2 zkv k,2 = 0, 



N-2 

-2Nz N - 1 a 2 + ( * k Vk,2 + N(N - l)z N ~ 2 ) ai + ^ z k V k , 3 = 0. 





N-2 N-2 
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We see that a\, a 2 are polynomials in z . By induction we get that any a s satisfies 
an equation of the form: 

N-2 

-sNz N - 1 a s + £ z^^iz- 1 ) = 0, 


where qk, s are already polynomials in z . This proves the lemma. □ 

Now we are ready to show that the operator A has coefficients Aj, which are 
polynomials in z~ x . Denote temporarily by r the degree of the polynomial 8, i. e. 
if 9{x) = 9 r z r + ..., then 9 r / 0. 

Lemma 2.4 The coefficients Aj of the operator A are polynomials in z . The 
degree of 8 r = m and 

A m = m , A m _i = 9 m -i, (2-7) 

Proof. Using the definition ([TJ) and applying the involution b to the relation 
9(x)K = KQ we get: 

AS = S9(d z ) 

As the coefficients of A are expressed as differential polynomials of the coefficients aj 
of S we get that Aj are also polynomials in z . Comparing the first two coefficients 



of the above equality we get also (2.7) . □ 



2.2 The string equation 

In this subsection we are going to show that for the bispectral operator L there 
exists another operator Q, for which the string equation ( |0.6D holds. This equation 
as well as other properties of the operator Q (with appropriate normalisation) would 
be crucial for our constructions. 

In what follows we would assume that the number m is divisible by N. This is 
not a restriction since we can always replace A by A^. We put m = Nl. 



Lemma 2.5 There is a natural number n such that: 

Q = K- l xd x nN+1 K, (2.8) 

is a differential operator. The operator Q is a solution to the string equation ( |0.6| ). 

Proof. Using the bispectral property one can write 

(&&L) m ~ l 6 = (-l) m - 1 b 1 ((adz N ) m ~ 1 A). 

Each application of the operator a,dz N to any differential operator P reduces its 
order by 1. Using the fact that the operator 

A = A m dr + A m _ 2 d 2 m - 2 + -, 
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where A m is a nonzero constant, we get that the operator 

(adz*)" 1 " 1 A = K m { S ,dz N ) m - 1 d™. 

is an operator of order 1. Now prescribing weights to z and to d z as follows: wt(z) = 
1, wt(<9 2 ) = — 1 we obtain that the right-hand side of the above identity has weight 
equal to (m — 1)N — m. This shows that the operator in the above equality has the 
form: 

(adz 7V ) m ~ 1 A = cz^-^-^dz + Cl z mN - m - N ,c^ 
In this way we get that 

Qi := (adL)" 1 - 1 ^ = 6i((-l) m - 1 (cz (m - 1)(7V - 1) a z + ciz mAr -" 1 - JV )) 

is a differential operator. Using the fact that m = Nl and that b\(z) = Ln we 
obtain 

_ ClL Ni-i-i) = cb^z^-^-^dz) = c6i(z n7V+1 a z ), 

where we have put n = l(N — 1) — 1. Now it is obvious that 

Q : = hiz^dz) = -((-lr^Qi - Cl L m - 1 - 1 ) 
c 

is a differential operator. The identity ( p.6[ ) is obtained by applying the bispectral 
involution to 

[z nN+1 d z ,z N ] =Nz N ( n+1 \ 

□ 

Corollary 2.6 For any positive integer i the following formula holds: 

(adLf(Q i )=i\N i L i( - n+1 \ (2.9) 
Proof. Assume that ( |2.9| ) is true for 1,2, ... ,i. Then 

(adL)^ 1 ^) =0 

Since ad L is a differentiation in the ring of differential operators with rational coef- 
ficients we can use the Leibnitz's rule: 

(adL) m (Q i+1 ) = (adL) m (Q\Q) = £ (* + ^ (ad L) i+1 ^ (Q*)(ad L)'(Q) 

The only nonzero term in the above sum is the one for j = 1, hence 

(adL) m (Q m ) = (i + ^.(adL)^)^ 1 . 
Now ( ^!s| ) follows by induction on i. □ 
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3 The infinite point 



The present section is divided into three subsections, corresponding to the basic 
results, which we shortly describe. In the first subsection we present the operator 
Q as a polynomial in L with coefficients - operators of lower order. The second 
result is a proof that the point oo is a regular singular point for the operator L 
(Theorem |0.3| ). In the last subsection we give an estimate of the degree n (of the 
string equation) in terms of the roots of the indicial equation at oo. All results will 
be used in performing Darboux transformations. Now we will fix the situation in 
which we are going to work. 

Definition 3.1 By O we denote the set of all functions that are holomorphic at oo. 
If we write a function V(x) from O as: 

11 

V(x) = x u {a$ + a\- + a 2 ^r H ), a / 0, z/ > 0. 

x ar 

then the number 

ord(y) := — v 

will be called order of V at oo and will be denoted by ord(V). 

We introduce also the ring 0[d] of all differential operators with coefficients from 
O. Obviously the bispectral operator L is from 0\d\. The only properties of L and 
Q relevant for our purposes in the present section are summed up in terms of the 
wave operator as follows: 

Definition 3.2 We say that the operator L £ 0[d] solves the string equation iff the 
following conditions are satisfied: 

1. There is a wave operator K = 1 + oiid^ 1 + • • • with coefficients from O for 
which LK = Kd N and ord(aj) < —1. 

2. There is an integer n > such that Q = Kxd£ N+1 K- 1 is a differential 
operator. We will call the pair (L,Q) a string pair. The minimal number n in 2. 
will be called the string number of L. 



3.1 Q as a polynomial in L 

For convenience denote by 1Z the differential extension of C[x] by adjoining the 



differential indeterminates yi,y2,- ■ - ,(see [Ka| for details). We endow the differential 
ring 1Z with graduation which will be useful in the sequel : for a monomial r = 
x no y^^ ■ ■ ■ y^ 8 ^ set wt(r) = uq — (n\ + i\) — • • • — (n s + i s ). This weight provides 
1Z with the structure of a Z-graded ring: 



n 



where lZ n is spanned over C by all monomials r € 7Z for which wt(r) = n. 

This graduation can be extended in a natural way to graduation of the ring of all 
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pseudo-differential operators with coefficients from the ring 7Z, by prescribing to the 
symbol of differentiation d weight wt(<9) = — 1. For convenience the last mentioned 
ring will be denoted by Psd7£. In this way a pseudo-differential operator 

P= ^ (l j i) J . (lj R 

j<m 

is homogeneous of weight n if for every j the coefficient a,j is a homogeneous element 
from 1Z n +j. If P is homogeneous then by wt(P) we will denote it's weight. We will 
need two lemmas. The proof of the first one being trivial will be omitted. 

Lemma 3.3 (i) Assume that P = d N + • • • is a homogenous pseudo-differential 
operator from Psd7£. Then P is invertible in Psd7£ and P^ 1 is homogeneous with 
weight —N. 

(ii) For every two homogeneous operators P\ and P2 from Psd 1Z with weights respec- 
tively n\ and ri2 their product P\.P% is also homogeneous and its weight is m + n%. 
□ 



Lemma 3.4 Assume that L y = d N + ■ ■ - and Q y are arbitrary homogeneous pseudo- 
differential operators from Psd 1Z. Then one can find an integer number n and 
homogeneous differential operators qo, q~i, . . . from lZ[d] of orders < N — 1 such that: 

Q y = q L n y + q x L n - x + • • • (3.1) 

More precisely, for any i = 0, 1, . . ., such that q~i ^ the weight of q~i is: wt(Q y ) — 
(n — i)wt(L y ). 

Proof. For given Q y we will show that q~Q and n can be determined uniquely and 
that they satisfy the properties stated in the lemma. After % is determined we move 
the term q~oL n to the left-hand side of ( p.l| ) and then in the same manner we can 
determine q\. Now it is clear that all % can be found successively and the lemma 
will be proved . 

Denote by m the order of Q v and divide m by N : m = nN + r, < r < N — 1. 



Multiply both sides of (3.1) by L n and compare the differential parts of the two 
pseudo-differential operators: 



qo = (QyL~ 



Combining this equality with Lemma |3.3| we obtain the statement of the lemma. □ 
Denote by PsdO the ring of all pseudo-differential operators with coefficients 
from O. To use the result of Lemma we need a ring homomorphism 

vr : Psdft — >Psd0 

defined as follows: take the unique differential homomorphism between 1Z and O that 
maps yi into aij, i = 1,2,... , where 014 are the coefficients of the wave operator K and 
then extend this homomorphism to homomorphism between Psd 1Z and PsdO by 



leaving d fixed. Now from the representation in Lemma 3.4 we can derive a similar 
one for the operators L and Q. 
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Lemma 3.5 Let L and Q form a string pair and n is the corresponding string 
number. Then one can find differential polynomials q~o, q%, . . . , q n from PsdTZ, such 
that if we set qi = Tr(q~i) then: 

Q = q L n + qiU 1 " 1 H \- q n . (3.2) 

The operators q~i are homogeneous. More precisely: % = xd x and if q~i 0, then 

wt(<7i) = — iN. 
The differential operator q n is not zero. 
Proof. Introduce the pseudo-differential operator 

K y = l + yi d~ l +■■■ € Psdft 

It is easy to check that L y = K y K~ l and Q y = K y xd^ n+1 K' 1 are homogeneous 
elements from PsdT?. with weights respectively: wt(L y ) = —N and wt(Q y ) = —nN. 
The definition of 1Z was given in such a way that ir{L y ) = L and 7r(Q y ) = Q. 
Applying Lemma 'iA with L y ,Q y we get: 



Qy — %L y + q\Ly 1 H , 

where each q~i is homogeneous with weight wt(gj) = wt(Qy) — (n + i)wt(L y ) = —iN. 
Map both sides of the last equality by it: 

Q = ir(q )L n + ■■■+ 7r(q n ) + Tr(q n+1 L y 1 +•••)■ 

Comparing the strictly pseudo-differential parts of the operators at the two sides of 
the above equality we see that: 

n(q n+1 Ly 1 H ) = 0. 

The inequality q n ^ holds because n was chosen to be the minimal number with 
the property that Kxd nN+l K~ 1 is a differential operator. □ 



3.2 x = oo is a regular singular point 



Here we give the proof of Theorem |0.3| , i.e. that the infinite point is regular for L. 
Take the smallest n for which Q = K xd^ n+1 K^ 1 is a differential operator. The 
idea is to assume that x = oo is irregular for L and then to assign weights to x and 
d x in such a way that the most irregular terms of L at oo have the highest weight. 
This weights will enable us to associate with each differential operator from 0[d] 
a (p, <r)-homogeneous polynomial in Y with coefficients Laurent polynomials in X. 
Following Dx and using ( p.6| ) we will get contradiction. 



We denote the ring of Laurent polynomials by C The definition of p and a is 
prompted by the theory of irregular points (see e.g. flBVfl )- Introduce the rational 
number: 

r = max(1 , 2 + 2^,.., 2 + ^) 
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(called principal level). It can be expressed as where n and r 2 are relatively 
prime. Well known fact is that x = oo is regular if and only if r = 1. Our assumption 
that oo is irregular point yields r > 1. The integers p = r 2 and a = r\—2r 2 represent 
the weights of x and d x respectively. They satisfy the inequality: 

p + a > 0. 

The next definitions are modifications of corresponding ones given by J. Dixmier 
[Dx|. In the first definition we endow the ring 0[d] (of differential operators with 
homomorphic at oo coefficients) with Z-graded structure. 

Definition 3.6 Assume that L = V d n + V^' 1 H \-V n is an arbitrary element 

of T>. For each term V(x)d x define its weight 

v P ,a(V(x)di) = p(ordV) + ai. 

Then the number 

v p ,a(L) := max Vp 1<T (V i d n ~ i ) 

0<i<n 

will be called (p, c)-order of L. 

The second definition associates to each differential operator from 0[d] a (p,cr)- 
homogeneous polynomial from £[Y]. 

Definition 3.7 Assume the notation of the previous definition and denote by I(L) 
the set {i G {0, 1, • • • , n}\v p ^(Vid n ~~' 1 ) = v p ^{L)}. The polynomial p G C\Y] defined 
as: 

p = ^a i X ordV >Y n - i , (3.3) 
iei 

where ai € C are uniquely determined from the expansion 
Vi = aiX ordVl + (lower order terms), 
will be called polynomial associated with L. 



The following two lemmas are also taken from [Dx]. Although the situation there 
is slightly different the proofs are essentially the same. We are going to prove only 
the first one. The second can be proven in a similar way. 

Lemma 3.8 Assume L\,L2 € 0[d] and p + a > 0. The polynomial associated 
to the product L1L2 is the product of the polynomials associated with L\ and L2 
respectively. The (p,o)-order of this operator is: v p ^{L\L2) = v Pj(T (Li) + v p ^(L<2). 



Proof. Set £ = d x . Then for the product of two differential operators we have: 

d k L 1 d k L 2 _ 

Qjrk Q x k 



fe=0 
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where : : is the normal ordering which always puts the differentiation on the right. 
Write L\ = ao^ Nl + • • • + aj\n > L 2 = ^oC^ 2 + ' ' ' + ^n 2 ■ From the definition of : : we 
have that 

: L X L 2 := £ a i b j ^ +N ^. 

0<i<N 1 ,0<j<N 2 

Each term in this sum satisfies the inequality v Pja {o.ibj^ Nl+N2 ~ t ^^) < v p ^(L\) + 
v p ,a(Li). The equality is possible only when i € I{L{) and j G I(L 2 ). On the other 
hand the coefficient in front of the highest degree of £ in: 

ie/(ii)j6/(ia) 

is a^ftij 7^ 0, where ii and Z2 are the minimal numbers from I{L\) and I(L 2 ) re- 
spectively. Thus this sum (which in fact is equal to the product of the (p,a) — 
polynomials associated with L\ and L 2 ) is not zero . The conclusion of this obser- 
vations is that the (p, cr)-polynomial associated with : L\L 2 : is the product of the 
polynomials associated with L\ and L 2 and also v P)CT (: L\L 2 :) = v Pja {L\) + v p ^(L 2 ). 
To finish the proof it is enough to use formula ( |3.4[ ) and the obvious fact that 
v Pt<T (d^Lx) < v Pj(T (L{) - ka and v p>a (d%L 2 ) < v p>a (L 2 ) - kp. □ 

Lemma 3.9 Consider again two operators L\,L 2 £ 0[d] and denote by fi , f 2 the 
polynomials associated with them and by n\ and n 2 their (p, o~)-orders. If the fraction 

"tttj- is not a constant and p + a > 0, then the polynomial associated with [L\,L 2 \ is: 
J 2 

dfidh dfidh ( , 

dYdX dXdY' [ ' ' 

For the (p, o~)-order we have a formula: v Pt<J ([Li, L 2 ]) = n\ + n 2 — p — a. □ 

In order to apply these lemmas to the string equation ( p.6[ ) we have to find the 
polynomials / and g associated with L and Q and their (p, <r)-orders v and w. This 
requires few auxiliary results, stated in the following two lemmas. 

Lemma 3.10 (i). The (p, a)-order of L is v = Na and the polynomial associated 
with L is: 

f = Y N + (at least one term). 

(ii). The (p,a)-order of Q is w = (nN + l)o~ + p and the polynomial associated with 
Q has the form: 

g = XYf n + ai (X, Y)f n ~ l + ■■■ + a n (X, Y). 

Proof, (i). Since v p ^ a (d N ) = Na the only thing we have to check is that ord(Vi)p + 
(N — i)a < Na for i = 2, 3, . . . , N and that equality is reached for at least one i. 
But this is obvious from the definition of p and a. 
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(ii). The polynomial g has the form: 



g(X, Y) = a (X, Y)f n + a x {X, Y)f n +■■■ + a n (X, Y), 



for some cij G Lemma |3.5| gives that a(,i = 1,2, ... ,n can have only negative 

degrees of X. But then the coefficient at the highest degree of Y in the polynomial 
g is not a constant, while the corresponding one in f w is 1. Thus the fraction 
is not a constant. Now from Lemma |3.9j the polynomial h associated with [L, Q] is: 

_ df_ dg_ _ df dg 
~dYdX dXdY' 

and v p>a (h) = v PjCr (f) + v Pja (g) — p — a = v + w — p — a. On the other hand the 
string equation Q0.6Q yields: v p>a (h) = (n+l)v. From the last two relation we derive 
the formula for w. To finish the proof it is enough to notice that v PiCT (qoL n ) = 
(Nn + l)a + p. a 

Lemma 3.11 Under the above notations g = XYf n . 



Proof. If h is the polynomial associated with [L, Q] then using Lemma 3.9 we have 
the following series of equalities: 

vi, v( d f d 9 d f d 9 , df dg 
= PX( WM -dxW ) = dY {W9 ~ aYdY9) ~ {vf ~ aY9Yf) 9Y = 

where we have used that for a (p, cr)-homogeneous polynomial / of (p, cr)-degree v 
the following identity holds: 



pXdxf + oYdyf = vf. 
Now the relation (|0.6| ) leads to: 

NpX.f n+1 = r w+1 g v+l dY (^j . (3.6) 

View / and g as elements in /CfK], where K, is an algebraic extension of the field of 
fractions of the ring C containing all the roots of the polynomials f(Y) and g(Y). 
Take a(X) to be a zero of / of order u > 1. Denote also by p the order of a(X) as 
a zero of g. Then comparing the orders of the terms at the both sides in formula 
( |3. 6|) we obtain 

vin + 1) + {W - X)V -(V+ l)p = OTdy. a{X) dy[-) . (3.7) 

We will treat the following 2 cases separately: 

easel. If wv ^ up,, then the right side of ( |3.7| ) is wv — vp — 1, hence p, = nv + 1. 
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case2. If wv = vp, then using the formulas for v and w we find p = —v = 

In both cases p > nu, which means that is a polynomial in Y. Now from Lemma 
0. 10 f n divides the polynomial a\f n ~ x + ■ ■ ■ + a n whose degree in Y does not exceed 



N-1 + N(n-1) < nN. But the degree of f n is exactly nN =>• ai/ n-1 H Va n = 0. 

□ 



Now we are ready to give the proof of Theorem 0.3 



Proof of Theorem D.3 . Put w = nv + p + a and g = XYf n in (|3.6[ ) and after 
simplifications we get a differential equation for /: 

Yd Y (f p+ n = (p+°)Nf p+ ' 7 . 

An immediate consequence of this equation is that / = c(X)Y N . But the choice of p 
and a was done in such a way that / = Y N + (at least one term) . This contradiction 
proves the theorem. □ 

The regularity of L imposes the following restrictions on the coefficients of the 
wave operator K. 

Corollary 3.12 Let L be an operator solving the string equation and K = 1 + 
a i d- 1 + ... is the wave operator defining the corresponding string pair. Then the 
order of the coefficient 014,1 = 1,2, .. . does not exceed the number —i, i.e. 

3.3 An estimate for n 

Here we want to estimate the number n from the string equation in terms of the 
roots of the indicial equation for L at 00. For us it would be convenient to write 
the indicial equation, using again the idea of the weights. But in order to have 



an analogue of Lemma 3.8 we have slightly to change the procedure of association 



polynomials to the elements of 0[d\. The next definition describes this process. 
Definition 3.13 Write every L G O[0] as 

L = V D N + ■■■ + V N -iD + Viv 
where D = xd x and assume also that: 

Vi = aiX Vi + (lower order terms) 

The number 

wt(L) := max ord(V;) 

will be called weight of L. The polynomial associated with L is from C[D] and is 
defined as follows: 



p(L) = Yl ai D 

i:ord(V;)=wt(L) 



N-i 
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In particular if the point x = oo is regular then p(X) = is explicitly the indicial 
equation (see ||). 

In terms of the above definition we can give the following corollary from Lemma 
U and Corollary [3~T2 . 



Corollary 3.14 Assume that (L,Q) is a string pair and n is the corresponding 
string number. Divide Q by L to derive 

Q = Q\L + q , 

where q is a differential operator of order not exceeding N — 1 . The weight of q 
satisfies the inequality: wt (q) < —nN. 

Obviously p(D)x l = x l p(D + i) for every polynomial p G C[D]. This observation is 
enough to make the following conclusion: 

Lemma 3.15 Assume that Li,L 2 £ T> are two arbitrary differential operators and 
denote by f\ and f% the polynomials associated with them. Then the polynomial 
associated with the product L\L 2 is: 

h{D + TNt{L 2 ))f 2 (D). 

The weight of the product is a sum of the weights of the two operators. □ 

Now we will assume that L is an operator solving the string equation. Denote by 
Ai, A2, • • • , Ajv the roots of the indicial equation of L at 00. The following very 
important fact, used in performing Darboux transformations, is the content of the 
next proposition. 

Proposition 3.16 Assume that L is a differential operator that solves the string 
equation. Then we can find numbers i and j such that: 



1 , X 

n < — I \i - A 



.1 



Proof. Let (L, Q) be a string pair. As in Corollary 3.14 divide Q by L 

Q = QiL + q 

Using that Q satisfies the string equation ( p.6[ ) we get 

Lq = LiL, (3.8) 

for some L\ £ 0[d]. 

Denote by /, g and h the polynomials associated with L, q and L\ respectively. For 
us the weight of q will be very important and will be denoted by w. 
Lemma 3.15 combined with ( |3.8| ) gives: 

f(D + w)g(D) = h(D-N)f(D). 

As a result we found that: /(Aj + w)g(Xi) = 0, for i = 1, 2, ...,N. Using the 
inequality: degg < iV — 1 one can find Aj for which g(\i) / 0, hence /(A« + w) = 0, 
i. e. Xj = X{+w, for some Xj. Applying Corollary |3.14| we get that w < —nN. This 
gives that 

nN <\ w |=| Xi — Xj I . 

□ 
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4 Z r - invariant ness of bispectral operators. 

Let Al be the ring of all differential operators commuting with L. We want to prove 
that if the rank of L is r then L is a Z r — invariant operator. The next lemma shows 
that it is enough to prove that A(z,d z ) is Z r -invariant. 

Lemma 4.1 If A is Zr-invariant then L is also Z r -invariant. 

Proof. It is enough to prove that the wave operator K is Z r -invariant. Assume that 
A is Z r - invariant. Then obviously = b{A) is also Z r -invariant. 

Now by induction on i we will see that the term otid~ l is Z r -invariant. Compare 
the coefficients in front of d~ 3 in the relation: 

0(1 + ai d~ l + •••) = (1 + "id" 1 + • • -)(©o + ©l^ 1 + •••)■ 

Comparing the coefficients in front of d° and d~ l one deduces that = Oo is 
Z r -invariant and that 0i = 0. Next assume that ot\d~ l , a2<9~ 2 , • • ■ , otid~ l are Zr- 
invariant and compare the coefficients in front of d~' l ~ 2 : 

«+2 / • o\ / • o\ 



The last formula together with the fact that is a Z r -invariant pseudo-differential 
operator and the inductive assumption give that q, + i9 _1_1 is Z r -invariant. □ 

The next lemma shows that the algebra Al consists of Z r -invariant polynomials. 



Lemma 4.2 Let L be an operator of rank r. Then Al is a subalgebra ofC[z 



r] 



Proof. Let P G Al- Put b\(P) = f(z) £ Al- From Lemma 1.9 we know that f(z) 
is a polynomial. Also the degree of f(z) is a number divisible by r. Assume that 
/ C[z r ] and also that the coefficient in front of the highest degree is 1. We can 
represent / as: 

/ = /o + /i, 

where /o S C[z r ] is formed from all terms of / whose degrees are divisible by r and 
fi = f — fo- The polynomial /o will be called the invariant part of / and /i the 
non-invariant part of /. Denote by no and n\ the degrees of fo and f\ respectively. 
Obviously no > n\ and n% is not divisible by r. The idea is to construct new 
polynomial / from A in such a way that the difference fiQ — h\ between the degrees 
of the invariant and the non-invariant part of / is smaller. After finitely many steps 
we will end up with a polynomial for which this difference is negative, which will be 
a contradiction. 

The polynomial / can be constructed as follows: let no = kr and N = pr set 
f ■= f p — z kN . Denote by /o and fi the invariant and the non- invariant parts of / 
and let no an d n\ be their degrees. Write the following chain of equalities: 

f = fP _ z N k = {fo + fi y_ 2 Nk = f P_ z kN + (p\ f p-l fi + ... 
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Since prig = kN and /q is a polynomial in z r we can conclude that fiQ < priQ — r. 
The above expansion together with no > n\ gives that h\ = n$(p — 1) + n\. Now we 
can prove that the new difference is smaller: 

no — n\ < pno — r — n\ = priQ — r — (p — l)no — n\ = n$ — n\ — r. 

□ 

Proof of Theorem CL4. It remains to prove the Z r -invariantness of A. Write A in 
the form: 



where 



r-l 



(4.1) 



i=0 



K{z\zd z ) = Y,Hi(z r )(z nN+l dz) n *- J 

3=0 

All Ajj are Laurent polynomials and is chosen in such a way that A^o 7^ 0, when 
Aj 7^ 0. We have to prove that all Aj, i = 1, 2, . . . , r — 1 are 0. Thus assume that at 
least one Aj 7^ 0. After applying the bispectral involution b\ on (4.1) we will get the 
following relation: 



no 



n r _i 



j=0 j=0 



r-l 
N 



(4.2) 



The idea is to construct an operator from Al whose image under the bispectral 
involution is not from C[z r ]. This will be contradiction with Lemma [4.2[ We split 
the construction of such an operator into two cases: 



easel, no < maxjni, ri2, . . . , n r _i}. 

Denote by p the maximal value of the numbers no, rai, . . . , n r _i and by / the set of 
all indeces i for which n, = p. Due to Lemma |2.l| 

(adLY(Q p ) = (p)\N p L^ n+1 \ 

hence one obtains the following relation: 

(ad Ly (9) = {p)\N p L^ n+1 ^ A ifl( L ^) L ^- ( 4 -3) 

iei 

Since the operator at the right-hand side commutes with L, it follows that the 
differential operator at the left-hand side is from Al- After applying the bispectral 
involution to ( |4.3| ) we get that: 

z P(n+l)NY^K, lfl (z r )z l 

iei 
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is an element from A^. This element is not polynomial in z r because the set I 
includes at least one index i £ {1, 2, . . . , r — 1}. 



case2. no > max{ni,n2, • • • , n r — i}. 
Now (|4.2|) can be written in the form: 



no — 1 rti 

^ - g n °Ao,o(i^) = E Q n °- J 'A 0)i (L^) + 2Q ni - j Aij(L*)^ + • • • + (4.4) 

j=l j=0 

n r -i 

j=0 



Applying (adL) n ° to the above equality we see (using Lemma 2.1) that it annihilates 
the operator at the right-hand side. Hence the operator 

(ad L) n °(Q n ° A<,,o(£*)) = (adLp(#) 



must be differential. Denote by Ni the number no(n + 1). Using again Lemma 24. 
i.e. that 

(adL)"°(Q no ) = noI^L" ^ 1 ) 



we see that after multiplying from the right both sides of (4.4) by L Nl the operator 



on the left-hand side will become differential. Denote this new operator by P. We 

ho 



re-denote A{jL Nl by Aj to avoid complicated notation. Thus the new relation has 



the form: 

■»»o-l/nno-7A_ .(T^r\ i Y^ n l 



P = q»o- JAOj . (jL1 v) + Q^A ld (LN) L N +...+ 

E]=o Q^-'K-i,^)^ ( 4 -5) 

We can repeat this procedure until no is reduced to a number smaller or equal to 
max {m, ri2, ■ ■ ■ , n r -i}- Then one proceeds as in case 1. □ 

5 Darboux transformations 

In this section we will gradually simplify the operator L by successive applications 
of Darboux transformations. Our goal is to obtain after a finite number of steps a 
Bessel operator. 



According to Theorem 0.2 the point x = oo is a regular singular point for the 



operator L. Assume also that L is a rank r differential operator. From Theorem 0.4 
we know that in this case L is Z r — invariant operator. Thus if we represent L as 

L = d N + V^- 1 + ■■■ + V N -id + V N 

then the coefficients Vi can be expanded as 

oo 

V i = -J2Vi, k x- rk (5.1) 
00 k=0 
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In what follows we need to split the set M = {Ai, A2, • • • , A^r} of roots of the indicial 
equation at 00 for L into subsets of equivalent modulo Z numbers. 

For an arbitrary set Mj denote by A the number in Mi with minimal real part. 
The next lemma is a version of a classical result (see e. g. B) and shows how one 
can pick an Z r -invariant function from KerL. 

Lemma 5.1 If X is the minimal number of a set Mi, then there is a function <f>\ 
from Ker L which can be expanded around 00 as : 

00 

A (s) = x x J2c k x- kr ,co = l (5.2) 

fc=0 

We omit the proof as it repeats the classical one. 

Given a function <p\ we construct a first order operator by setting 

P _« fix 

P X — Ox - — ■ 

(f>X 

Then the operator L can be factorised as L = Q\P\ and after we perform the 
Darboux transformation 

L = QxPx — * L = P X Q X (5.3) 
the new operator L will have the following properties: 

Proposition 5.2 Assume that the operator L solves the string equation with a Zr- 
invariant wave operator K = 1 + ai<9 -1 H , ord(aj) < —i. Then 

(i) every operator which is obtained by a Darboux transformation described above 
also solves the string equation; 

(ii) if{Xi, A2, . . . , Ajv| are the roots of the indicial equation at 00 of L and A = Aj 
is the number with minimal real part from some Mj then the roots of the indicial 
equation at 00 of L are A^ = A& — 1 for ,k^io and A.; = Aj + (N — 1). 

Proof. Put 

K = P x Kd-\ (5.4) 
Now we will check that L, K also satisfy the conditions of the lemma. Let's check 



the first condition of Definition 3.2 



LK = PxQxPxKd- 1 = P\LKd~ l = P x Kd N d~ l = Kd N . 

Further denote by Q = PxQQx and note that the following sequence of equalities 
holds: 

QK = PxQQxPxKd- 1 = PxQLKd- 1 . 
Using the equalities LK = Kd N and QK = Kxd nN+l the last relations give 

qk = p x Kxd nN+l d N - 1 = Kdxd N ^ = kxd^ N+l + kd^ N . 
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Now it is clear that Kxd^ n+1 ^ N+1 K 1 = Q — L n+l is a differential operator. 

Denote by g the polynomial associated with P\ and by h the one associated with 



Q\. Obviously g(D) = D — A and it has weight —1. Then using Lemma 3.15 we get: 

h(D - l)g(D) = (D - X 1 )(D - A 2 ) ...(£> - A,v) 
g ( D — (TV — = (D — Ai)(D - A 2 ) ...(£> - Ajv) 

^,From these equalities we get the second assertion in the lemma. □ 
After this proposition we are close to our final goal. 

Proposition 5.3 Let L be a bispectral operator with coefficients satisfying ( p.4[ ) . 
Then by finitely many Z r — invariant Darboux transformations we can transform it 
into a Bessel operator. 

Proof. We will perform Darboux transformations in the following way: start with 
Lq = L. Choose an index i, if there is any, for which the difference between numbers 
in Mi with maximal real part and with minimal real part exceeds N. Denote by 
A the number in M, with the minimal real part, set Pi = P\ and factorise L as 
L = R\P\ then the Darboux transformation will be 

Lq = L = R\P\ — ► L\ := P\R\ 

According to Proposition [|[^ the sets M® := Mj will be transformed into sets Mj 
for which the difference between the numbers with maximal and minimal real parts 
are the same for % ^ j. When i = j there are two cases: 

easel. There is exactly one number in Mi with minimal real part. The differences 
between the numbers in Mj are integer. Thus there is a well defined ordering: A > \i 
iff A — n > 0, in fact A — fj, = ReA — Re^u. Having in mind this remark the elements 
of Mi can be ordered as 

A < /xi < ... < fj, s . 
Now the assumption about Mj means 

A+(AT-1) <n s -l. 



Due to Proposition 5.2 in the new set Ml the following inequalities must hold: 
minM^ 1 > A, maxM? = fj, s — 1. Hence, the difference between the maximal and 
the minimal number is reduced at least by 1. 



case2. In Mj there is at least two numbers with minimal real part. Then the above 
Darboux transformation decreases the number of the roots with minimal real part 
at least by 1. 

After finitely many Darboux transformations we obtain an operator L m such that if 
M m is the set of roots of the indicial equation at oo and Mj 71 are the corresponding 
subsets modulo Z for M m , then 

max Mf - min Mf < N. (5.5) 
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But again from Proposition 5.2 it follows that there is an operator 

K m = l + aid" 1 + ■■■, 
such that L m K m = K m d and there is an integer n > for which 

Q m = K m xd nN + x K^ (5.6) 

is differential. The minimal n with this property, according to Proposition 3.16| 
satisfies the inequality: 

n < — (max M. m - min Mf ) . 
- N y 3 3 ' 

Using (5.5) we see that n must be zero. Put in Q5.6D n = and compare the 
differential parts of the operators at both sides to conclude that 

Qm — xd x - 

Now comparing the coefficients at d 3 at both sides of the string equation Q0.6| ) with 
n = we obtain the equation 

-xVj + jVj = NVj. 

Integrating it, we obtain that 

Vj = VjX~ N+j , Vj eC. 
This shows that L m is a Bessel operator. □ 

Proof of Theorem |0.2| . It remains to show that the chaine of the above Darboux 
transformations can be replaced by one monomial. First we represent the chain by 
following graph: 

Lq = R\P\ — > L\ = P\R\ = R2P2 —> L2 = P2R2 = R3P3 L m = P m R m . 

If we set A = R1R2 ■ ■ ■ R m and B = P m P m _i ... Pi then obviously: 

L m = = R 1 P l R 1 P 1 . . . RxPt = R 1 L r [ l - 1 P 1 = AB 
and for the Bessel operator Lri := L m 

Lgi = P m R m P m R m . . . P m R m = PmJ j ra —\Rm = BA (5-7) 

The Darboux transformations do not change the rank of the operator. Thus the rank 
of Lri is r. If r < N then according to Lemma 1.6 there is a monomial Darboux 
transformation which transforms Lri into Lr, where Lp is some Bessel operator 
of order r and rank r. But the monomial Darboux transformations connecting 
Bessel operators are transitive. Thus there is a monomial Darboux transformation 
connecting L and Lp. The only thing that we have to prove is that the operators 
A and B from (5.7) have rational coefficients. To prove this we need the following 
lemma. 
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Lemma 5.4 Assume that P G 0[d] is an operator with holomorphic at oo coeffi- 
cients. If P divides from the right some power Li of a Bessel operator 

L p = x- N (D-(3 1 )...(D-(3 N ), D = xd x , 

then the coefficients of P are rational. 

Proof. Let n be the order of P and 

7 = p d = (ft, A + N, . . . , ft + (d - 1)N, . . . , /3 N , I3 N + N, . . . , (3 N + (d - l)N). 

First we prove that Ker P has a basis of elements /j, i = 1, 2, . . . , n of the form: 

n 

fi = x» Y^PiA^xy ^in + 0, (5.8) 
j=0 

where are polynomials. In general every / € Ker P can be written as: 

s 

f = T,fi (5-9) 

with /j having of the form given by ( |5,8D and 7, —7,- Z for i ^ j. The analytical 
continuation around the infinite point defines the monodromy map: 

Moo : Ker P — ► Ker P. 

If an element / = J2i=i fi as i n ( |5.9| ) and J5.8|) is in Ker P, then 

s n 

i=X j=0 

is also in Ker P. 

Let s be the minimal number for which there is an element / as in (|5.9|), where 
none of the terms fi is in Ker P. From all such functions from Ker P with minimal 
s take one for which the number: 

minjrj | i = 1, 2, . . . , s} 

is minimal. We can assume that r s = min{rj | i = 1, 2, . . . , s}. Then in the following 
element from KerP: 

s n 

f _ exp(-27r7 sV / ^l)M 0O (/) = £ /; = x^ ^p-,(x)(lnx)^ 

i=l j=o 

either the term f s vanishes (when r s = 0) or the number f s = r s — 1 is less than r s . 
In both cases this is a contradiction with the choice of /. 



28 



Having in mind the basis from ( |5.8| ) the action of the operator P can be written 
as (see 

p , Wr(fi,f 2 ,...,fn,(f>) /, 

P( f >= mt it r T~r- ( 5 - 10 ) 

Wr(f 1 ,f 2 , ...,/„) 

Note that each derivative /-^ has the form = x" li Fik{x, lnx), where Fi^(X, Y) E 
£[y] is a polynomial in Y with coefficients - Laurent polynomials in X. Hence formula 
( |5.10| ) gives: 



x 



71+72H h7s 



x~n+i2+-+is p n (xMx), 

where Fi E £[Y]. Thus the coefficient Cj in front of d l is 

Fi(x, lnx) 
F n (x,lnx)' 

Since q S O the monodromy map preserves q. Hence 

F i (x,lnx + 27r v /r l/) _ Fi(x,\nx) 
F n (x, lnx + 2■K^J^\l) F n (x,lnx) 

for every integer I and also for every / G C since the above equality is equivalent to an 
equality between polynomials. Using that x and lnx are algebraically independent 
we get: 

Fj(X,Y + l) _ Fj(X,Y) 
F n {X,Y + l) F n (X,Y) 

which on the other hand is equivalent to 

Fj(X,Y + l) = F n (X,Y + l) 
Fi(X,Y) F n (X,Y) ' 

Take the derivative with respect to I and set I = 0. Then one sees that Fi(X,Y) = 
c(X)F n (X, Y). After putting first Y = lnx, X = x and then Y = 0, X = x it follows 
that 

Fi(x,\nx) Fi(x,Q) 

Cj(X = — — ; : = ClX = — — r 

iV ' F n (x,lnx) V ' F n (x,0) 
is a rational function. □ 



6 Proof of the characterisation theorem. 



Essentially the proof of Theorem D.5 has already been performed in the previous 
sections, as well as in [BHY1, BHY5]. Below we sketch a plan how to pick the pieces 
of the proof from these sources. 
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Proof of Theorem p.5| . The implication 1) — > 3) is the content of Theorem 0.2. 
Next we consider 3) — > 2). Here we use the Definition L4 for monomial Darboux 
transformations. If Lp is a Bessel operator then one factorises Ua as 

= QP, (6.1) 

where the operator P acts on tp in the following way: 



P 



Wr{h,f 2 ,...,f n ^) 
Wr{h,f 2 ,...J n ) 



(6.2) 



and the functions f\, . . . , f n have the structure prescribed in Definition 1.4. Having 
in mind the type of the kernel it is obvious that the operator P has only regular 
singularities. But then the same is true for the operator Q whose coefficients are 
computed by induction from the fl6.1|). T hen the same is true for the product PQ. 
At the end by the main result in |BHY1| 1 the latter operator is bispectral. 

The implication 2) —* 1) is trivial. The equivalence of 3) and 4) is the content 
of [|BHY5|1 . We briefly describe it. 

First, we recall the definition of Wi +oc , its subalgebras Wi +00 (N) and their 
bosonic representations introduced in [BHY1]. The algebra Woq of the additional 
symmetries of the KP-hierarchy is isomorphic to the Lie algebra of regular polyno- 
mial differential operators on the circle 

V = span{z Q df | a,P E Z, (3 > 0}. 

Its unique central extension |KR| will be denoted by Wi +oa . This algebra gives the 



action of the additional symmetries on tau- functions (see [ASvM, |OS| ). Denote by 
c the central element of Wi+oo and by W(.A) the image of A £ T> under the natural 
embedding V •— > W\ +00 (as vector spaces). The algebra Wi+oo has a basis 



c 4 = w(- 

The commutation relations of W\. 
generating series [KR 



■z 



l+kal 



Uez, i > o. 



can be written most conveniently in terms of 



W(z k e 



k „xD, 



),W{z m e yD * 



(e J 



)W(z 



k+m e (x+y)D z 



+ $k 



e xm _ e yk 



1 



C, 



(6.3) 



where D z = zd z . 

^From the theory of KP-hierarchy it is well known that each operator L or its 



wave function (1.2) defines or can be defined by the so called tau-function, which is 
a function r(ti, . . . ,t n , . . .) in infinite number of variables t n ,n = 1, . . .. We denote 



the tau- functions of the Bessel operators Lp by Tp. In |BHY2] a family of highest 
weight modules Aip over Wi +OQ has been constructed, using as a highest weight 
vector rg. We briefly describe them. 

Introduce the subalgebra Wi +ao (N) of Wi+oo spanned by c and J l kN , l,k G Z, 



I > 0. It is a simple fact that Wi +00 (N) is isomorphic to W\ +00 (see 
put 



Mp = span{j l k \ N ■ ■ ■ J l k p pN rp 



h < ... < k p < o}. 



|). Now 
(6.4) 



The main result of [BHY5] can be summed up as: 
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Theorem 6.1 If an element in a module Mp is a tau-function then the correspond- 
ing operator L is a monomial Darboux transformation of some Bessel operator Lw 
(with eventually different (3' ). If an operator L is a monomial Darboux transfor- 
mation of a Bessel operator L@ then the corresponding tau-function belongs to the 
module M.p. 

Obviously the above cited theorem gives the equivalence between 3) and 4). □ 
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